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ABSTRACT 

The power spectral density of a pulse-position modulated (PPM) 
waveform is determined in part by certain statistics of the modulating 
process. In this paper a mathematical model is developed to determine 
the effects of three different modulating processes upon the power spec- 
tral density. The processes considered are the Gaussian, the exponen- 
tial, and the uniform processes, with statistical independence assumed 
between the sample values that modulate the PPM waveform. The 
Gaussian process is also developed for the case in which statistical in- 
dependence may not be assumed; that is, a correlation exists between 
the sample values of the modulating process. The model will also 
determine the power spectral density of pulse-amplitude modulated 
(PAM) waveforms, or for a combination of PAM and PPM, 

The mathematical model is reduced to a computer model and 
computed results are displayed. The design of a pulse-position modu- 
lator is shown, and experimental measurements of power spectra are 
compared with the theoretically determined spectra. The theoretical 
determination of power spectral density and the actual measurement of 
spectra of finite-time samples are discussed. 

The major portion of the research reported on was conducted by 
the writer during his Industrial Tour at Sylvania's Electronic Defense 
Laboratories (EDL) in Mountain View, California. 

The writer gratefully acknowledges the helpful as sistance and 
encouragement received from both his thesis advisor, Professor George 


H. Marmont, and from Dr. Paul O. Scheibe, formerly of EDL. 
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hs Introduction 

The objective of this report is to relate a theoretical approach 
for the determination of power spectral density, to show computations of 
power spectral densities performed for a particular class of signals, and 
to compare experimental measurements of the power spectral densities 
to that predicted by theory. 

The power spectral density of a pulse-position modulated (PPM) 
waveform typically consists of a continuous component and discrete fre- 
quency components, The relative amount of the total power contained in 
each of the component types is determined by the modulation process, 
The distribution of power in the spectrum is determined by the modula- 
tion process, the pulse shape, and the average pulse repetition rate. 

In this study, a mathematical model is developed which determines 
the power spectraldensity of PPM waveforms. The mathematical model 
is then reduced to a digital computer program. The computer program 
permits the computation of the power spectral density of a PPM wave- 
form that is modulated by a uniform process, a Gaussian process, or an 
exponential process. The modulation process that results ina statistical 
correlation between pulses within the pulse train with respect to their 
displacements (in time) is developed for the Gaussian process only. Sta- 
tistical independence between pulses is assumed for the other processes 
considered. 

The pulse shape is assumed to be rectangular for the computer 


model. The pulse width and the average pulse repetition rate are variable 


parameters of the program, The pulse amplitude is also a variable param- 
eter, thereby allowing computations to be performed for combined pulse- 
position and pulse-amplitude modulated waveforms. By setting the pulse 
position modulation equal to zero, the computation of the power spectral 
density of a pulse waveform that is only pulse-amplitude modulated may 
also be accomplished with the computer model. 

A discussion of the use of finite-time interval measurements, 
using a periodogram technique to give estimates of the power spectral 
density of a process, is included inthis report. The periodogram tech- 
nique is used to determine the power spectral density of a PPM waveform, 
Since the characteristics of the experimental PPM waveform are also 
used as an input to the computer model, a comparison of the theoretical 


and experimental results is included. 











rae A Discussion of Spectral Analysis (General) 

The purpose of a communication system is to convey information 
in some form from one place to another, The system would normally 
include a transmitter, a transmission medium, anda receiver which 
recreates a reasonable representation of the original message. The 
message creating the transmitted signal will most likely be one ofa 
class of stochastic phenomena; that is, the message is a random func-~- 
tion of time. Also, the noise accompanying the signal (introduced by 
the communication system) will normally be a stationary Gaussian ran- 
dom process. The analysis of such signals is difficult unless one re- 
sorts to the use of the statistical properties of the signals. The power 
spectral density furnishes us such an analysis technique, as it is rela- 
ted to the statistical parameters of the signal. In particular, the power 
spectral density is related to the covariance (or autocorrelation) function. 

In the development of spectral analysis techniques, one normally 
considers the random signal to be a stationary process, A stationary 
process is one in which the statistical parameters of the signal are the 
same everywhere, independent of which long interval over which they 
are measured, A similar concept is to consider a process as an en~ 
semble of signals (an ensemble is a set of similar things), If the statis- 
tical averages of the ensemble measured at one instant of time are the 
same for any other instant of time, then the process is also said to be 
stationary in time. 


If a single process of the ensemble of signals has statistics over 





a long period of time that are the same as the statistics of the ensemble 
at a given time, then the process is said to be ''ergodic". An ergodic 
process will always be stationary, but a stationary process does not 
necessarily imply ergodicity [6]. 

In this paper all processes will be assumed stationary unless 
otherwise stated. 

To arrive at the statistical parameter used to obtain the power 
spectral density, consider the random ergodic process y(t). <A typical 
member of the ensemble is y(t). The ensemble averaged correlation 
function (covariance function K(T))is defined as the second joint moment, 


or 


i) = EL eent ? 


] 
where E denotes the expected value taken over the ensemble. The auto- 


- : n . : 
correlation function R (T) of a typical member of the ensemble is 


lim I 


T 
R' (tT) _ T-+o DT do y (t)y (t + T) dt 


In general, 
K(t) = E{R(1)}, 
where the E{R"(1)} here is the autocorrelation function averaged over the 


The statistical definition of covariance would require that the 
square of the mean (uy“) be subtracted from the above expression, It was 
decided, however, to include the mean in this investigation, thus permit- 
ting a simplified development that will handle signals that possess a non~ 
zero mean. 


4 





entire (infinite) ensemble. This average should normally be performed 
before the time average of R(T), 


However, for the ergodic case 
Te KU (probability 1). 


The mathematics relating the power spectral density to the autocorrela- 
tion (or covariance) function were first developed by Wiener and 
Khintchine. The Wiener-Khintchine Theorem states that the power 
spectral density S(f) is the Fourier transform of the covariance func-~- 
tien, OT 

OO : 

S(f) = 2 | K(t)e ar, f 20 

~0O 
Note that in the above expression S(f) is defined for positive frequencies 
only, hence the factor of two in the equation. This is obviously the more 
practical case, and is permissable since S(f) is an even function (as is 
R(T)). One can obtain the autocorrelation function by taking the inverse 
Fourier transform of the power spectral density (again assuming ergo- 


dicity), Thus 


Jeg) 


Ne 


OO 
| S(f)cos wtdf (probability 1) 
O 
and 
OO 
S(f) = 4| R(T)cos wTdT (probability 1) 
O 


It should be mentioned that while a signal has a unique power spectral 


ae) oan 





__— -_ 








density, the converse is not true, All phase information in the signal 
is lost in S(f), thus the original signal could not in general be recov~ 

ered from its power spectral density. Furthermore, a particular S(f) 
theoretically represents the power spectral density for a multitude of 


Signals. 








3% Finite-Time Interval Spectral Analysis 

The developments of the previous section followed the typical 
theoretical approach to signal analysis in that infinite time intervals 
were assumed. In practical analysis situations, such intervals are of 
course always finite. For some processes, a finite investigation in- 
terval can be considered adequate if the interval is large with respect 
to the fluctuation time of the process under investigation. Since finite 
time measurements of spectra of PPM waveforms will be compared 
with the results from the mathematical model, a discussion of possible 
sources of error due to analysis of the finite-time interval observation 
of the signal is included. 

The problem of finite time interval analysis of a stationary ran-~ 
dom process can be stated quite simply as follows. A sample function 
of the random process y(t) is observed during an interval t = -T/2 to 
t = T/2. Will this observation contain enough of the characteristics of 
the process to give a reasonable estimate of the power spectral density 
of the process? Note that the random process under consideration is 
stationary. The stationarity of the signal becomes more meaningful 
(and necessary) when considering finite-time analysis. 

To develop the theory for finite-time analysis, again consider 
the ensemble y(t), except now let y(t) vanish outside the interval (-T/2, 
T/2). An amplitude spectrum for a member of the ensemble is 


n n jWt Tre n jWt 
Y"(£).. = i y' (t)pe dt = | y (tle > dt 
~CO -T/2 








Now the finite-time power spectral density of the ensemble is 


sity = B{Z O_O } = EE LI g Pf. 


* , 
where denotes the complex conjugate. The true power spectral den- 


sity S(f) is defined in terms of the ensemble average over an infinite 


time interval (and is equivalent to the definition on page 5). Specifically, 


lim 
i co 


_lim 2 


Sif) = Sif), = E{|y"0.1°} 


dh o'e)eih 


where the ensemble average is carried out before the limit is taken. 
Again, consider the sample function of a single random process. 


The finite time transform is equivalent to Y'(f)_, of the preceding para- 


T 


graph, and will now be uSed as an approximation to the Fourier trans- 
form of a single stationary process, A finite power spectral density of 


the single random signal may now be defined. 


n Z n 2 
Sra) ee | ate 


2; 


S(f, T) is called a 'periodogram", and has received considerable 


attention due to the fact that it is one type of practical measurement 


that can be made on real processes (signals), Observe that the power 


] 
Middleton, D. An Introduction to Statistical Communication 
Theory. McGraw-Hill, 1960, 143 


Abramson, N. and Farison, J. On Statistical Communication 
Theory. Technical Report No. 2005-1, SEL-62-078, Stanford Univer- 
sity, August 1962, 30-33 


3 
Grenander, U. and Rosenblatt, M. Statistical Analysis of 
Stationary Time Series. John Wiley & Sons, 1957, 91-94 
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spectral density for the time limited ensemble is the expected value of 


the periodogram over the entire ensemble, or 
n 
S(f),, = E{S™(£,T)} 


and 


im 


Se all E{s"(£,T) } 


Thus, one measure Of the reliability of the periodogram is the fact that 
as Too, the mean value of the periodogram does equal the power 
spectral density. However, as shown by Abramson and Farison, for 
some very common random processes, the variance of S(f, T) does not 
go to zero as Too, : The fact that the periodogram (being a random 
variable) has the desired mean value but may vary extensively about the 
mean for certain processes has caused some skepticism about its use 
as a spectral density estimator [5]. 

It has been shown that the periodogram is not a consistent esti- 
mate of the power spectral density [10]. However, the periodogram 
may be used to give consistent estimates of the power spectral density 
at specific frequencies f.. This is accomplished by the use of weight- 
ing functions, or "spectral windows" about the frequency f.. A possible 
method then of obtaining a consistent estimator of the power spectral 
density would be to have many adjacent windows, and at the end of the 
analysis time T observe the response of all windows simultaneously. 

A smoothing between the point estimates observed would naturally be 


desirable, 





Periodogram calculators fall into three general classes; that is, 
the filtering, computing and smoothing may be performed in three dis- 


tinctly different ways. These are 


le Parallel filter banks, 
oe Sweeping analyzers, 
33 The coherent memory delay integrator. 


The parallel filter bank technique is, as the name implies, a bank 
of narrow-band filters, each with a common input (the signal to be ana- 
lyzed). At the output of each filter is an energy detector, which is 
periodically sensed by a sweep eee that is connected to the vertical 
deflection plates of an oscilloscope. The parallel filter bank technique 


is depicted in block diagram form in Figure 3-1. 


Input Spectral window ] Energy 
tO ft Detector Sweeping 


Readout 
Device 





| 
| 
| 
! 
| Output 
| 

| 


| 
| 
| 
| 
| 
| 
| 
| 


Spectral window Energy 
for fy Detector 
Figure 3-1, Parallel Filter Bank Periodogram Calculator. 


The number of filters (N) may vary with the frequency range (analyzer 


bandwidth) and the resolution required. The bandwidth of a filter is 


10 





approximately the inverse of the rise time. Obviously, the bandwidth 
of the filters determines the resolution of the device. The "analysis 
time", or Period T, using this method is approximately the rise time 
of the filter, From this, one may observe that the resolution and the 
analysis time are essentially the inverse of each other. Consequently, 
if one designs an analyzer for fine resolution, then the signal must be 
observed for a longer period, and the converse is also true. This im~ 
portant relationship holds for periodogram calculators in general. 

The sweeping analyzer is the more common general laboratory 
type of device. With this device, the input signal is mixed with a sweep~ 
ing local oscillator, the output of which is filtered through a single 
"spectral window''. The signal is then detected and displayed on an os- 
cilloscope. A block diagram for the sweeping analyzer is shown in 


Figure 3-2. 






Narrow 


Band 
Filter 








Input 





Detector Output 


Sweepin 
Local = 
Oscillator 


Figure-3-2. Sweeping Analyzer Periodogram Calculator. 


A disadvantage of this type of analyzer is the adverse effect of the tran~ 


sient condition at the input to the filter while the local oscillator is 


BE 





sweeping. If the sweep rate is too rapid, the filter may continue to ring 
due to this transient effect, even though the signal is no longer present 
at the input to the filter, and the resolution is no longer the bandwidth 
of the filter, 

The third type of periodogram calculator to be discussed is the 
coherent memory delay integrator. This device is somewhat sensitive 
in its operation, but the concept is unique and merits consideration, 


The block diagram of the device is shown in Figure 3-3, 


Readout after 


Input Delay 
T 


Unity Gain Storage Loop 


Am a 
(Phase 





Shifter) 


Phase Advance Signal 
: an), 
cos | 
Figure 3-3, Coherent Memory Delay integrator 


The local oscillator used for the phase advance signal must be precisely 
matched to the delay of the delay line. Furthermore, the loop gain must 
be unity. With these optimum conditions assumed in being, the response 
of the analyzer to a single frequency is a main lobe, with side lobe pat- 


terns very much like a (sin ube curve. The side lobe patterns may 


i2 


be removed in some instances by a process of selective "weighting", 
in which the input signal is time weighted during the analysis interval. 
As an example, a cosine weighting function is accomplished by ampli- 
tude modulating the input signal. 

The output is produced by a linear addition of a large number of 
sinusoids differing from each other in frequency by small increments. 
Phase coincidence of all signals present in the storage loop occurs once 
(for a given frequency) during the observation interval. The time of 
this phase coincidence is linearly related to the input frequency, thus a 
readout of the response after the proper number of loop recirculations 
yields the desired periodogram, The one operational analyzer of this 
type with which this writer is familiar is known under the trade name 


"Simoramic'"', 


V3 








4. The PPM Spectral Model 

Pulse-position modulation (PPM) is a form of pulse time modu- 
lation in which the modulating process causes the position (in time) of 
the pulse to vary relative to its unmodulated time of occurrence, The 
relation of the pulse to the modulation process is normally determined 
by the sampling technique. The two more common types of sampling 
are natural sampling and "uniform" sampling [2]. For natural sam- 
pling, the occurrence of the pulse is determined by the instant at which 
the leading edge of a sawtooth sampling wave and the modulating (mes-~ 


sage) voltage are equal. 


yo :~S:*C«sMe Sage waveform x(t) 







sawtooth sampling wave 


== f 


resultant PPM 


' = waveiorm yi(t) 
tl {1 fT) fa 


Figure 4-1, Generation of PPM by Natural Sampling. 


When uniform sampling is used, the modulating waveform is 
sampled at equally spaced time intervals, Each pulse in the PPM wave- 
form is then delayed in time from its unmodulated position by a time e 


which is proportional to the sampled voltage. 


14 





4 message waveitorm x(t) 


sampled value x, 


resultant PPM 


Figure 4-2. Generation of PPM by Uniform Sampling. 


If one desires to relate a specific modulating process to the re- 
sultant power spectral density of the PPM waveform, it is necessary 
that the pulse positioning process € be expressable in terms of the mod- 
ulating process x, This was done for the case of uniform sampling; 
however, the attempt to do this for the natural sampling case proved 
fruitless, The resulting expression for the natural sampling case was 
in the form of an unwieldy series of expansion which could not be reduced 
to a closed form nor was of any value in its series form. Therefore, 
the spectral model development will be based on a single channel PPM 
system using uniform sampling. 

The zero order hold is employed by the uniform sampling tech- 
nique just described. The model can be adapted for use with other sam- 
pling schemes, For a Gaussian process, one can use a higher order 
hold in the model, ; as the model will be developed to handle the case .~ 

An example would be a technique which requires that the value 
a be linearly related to the average of the previous N samples of 


5 





for which ae 2 are correlated Gaussian bivariate random variables 
for a finite m-n. The correlation for this hypothetical case can be in-.: 
troduced by the sampling technique, or it may be directly related to the 
correlation of the sampled values of the modulating process, or both. 
Let the modulating waveform x(t) vary from zero to A volts. 
Define the maximum allowable deviation of the pulse from its unmodula-~ 


ted position as T/L seconds, The time shift é, is then 


as T 
emis 


O 
a Wes 


The pulse time shift €, is thus linearly related to the sampled amplitude 


k 


of the message waveform x(KT.). The variance of ¢, is also linearly 


k 
related to the variance of the message waveform samples x» Or 


: 1% iL Ze 
“ E{e a E{—° x = } fe) 
€ Hef 7 ee in ee 


2 (Tee 2 
e- \LA <—— 


The standard deviation (or RMS deviation from the mean) of the pulse 


Q 
tl 


time shift ¢, is 


A more meaningful expression in determining what the actual 


RMS deviation is relative to a normalized interval will be defined as a 





A general expression for the power spectral density of a PPM 
waveform will be shown, The process x(t) producing the modulation 
will be considered stationary. The development of the equations will 
follow very closely that of Middleton [8]. 

A pulse function U(t ~ nT. - €) shall be defined as having any 
desired shape (with maximum value of one) and existing only for a dura~ 
tion u beginning at time t = nT’. ts c. Outside of this specified inter~- 
val, the pulse function is zero. Define y as the peak amplitude of the 
7 member of the ensemble of PPM waveforms in the time interval 
nT 1 Sus <t< nT. at; a +u. A member of the truncated ensemble 


Yqylt) may now be written: 
N 


1 ser 
Yn (t) = )) ie Uit - nT) ~ € 
n==-N 


i 
) 


n 


The autocorrelation function of the PPM waveform 


i igen i 
is R(T). 
2 
4 lim 1 i i 
Ret) = 2 | ye tye at 
5 
co ak _T/2 N N 


ihe 1 : 7 | 
~ Noo (2N + 1)T. Ym Yn 


m, n=-N 
rs i i 
| U(it -mT -e )U(t+ T=nT +c jdt 
ann O mM O 1G 


Now let Xy(f) be the Fourier transform of U(t). 
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x (f) = | U(t)e ) dt 


-OO 
or 
sad }wt 1}wWO 
| U(t -a)je ? dt = X (fe ley. 
=OCO 
Also 
2 ‘wt 
U(t) - | Xy(£)e” df 
-0O 
or 


CO ‘ 
U(t - B) = | x, (fer Pag 
-OO 


Letting a = nT. + co - T and $8 = mT. + cn in the above equations, we 


have 
N CO 
it lim ] )) i 1 
i ee — a a ~ 
SN, N-oo (2N + 1)T Yn Ym | Ut - @) 
Oom,n=-N =00 


CO ‘ 
| x, (fe asat 


-OO 
N 
_ |lim 1 ) i (x (f) 67 508 
~ |N-00 (2N + 1)T Yn Ym U 
om, n=-N aoe 


= jwt 
| U(t - ae” dtdf 


-00 
N 
- ie — ) i i x (£97508 
~ IN-oo (2N + 1)T Ym n U 
L om, n=-N a 
*  Jwa 
X_,(f) e° df 
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m,n=-N 


Now is it necessary to determine the covariance function from the above. 


Recall that the covariance function is defined as 
K(T) = E{R(7) } 


Additionally, recall that it is necessary to perform the ensemble aver- 


age prior to taking the time average of the autocorrelation function. 


Thus 
CO 
1 | 2 -jwTiliee: l 
T= — ——~—__ 
ESCO) | IX, (8) | . Noo (2N + 1) 
O ~=0O 
N .  . jwie v Le = jw(n-m)T 
E{ ¥ a 1 n m : °f df 
Uae vs 
m,n=-N 


Now observe three important facts. The expected value of a finite ser- 





ies 1s being taken, such that an vanish for lm |, \n | 2N. Since 


the process is stationary, the statistical relationship of any two random 
variables displaced the same distance kT in time is the same for all 
such pairs of random variables. Finally, there are 2N + 1 such like 


processes in the double series. Define k =m -n inthe series, and the 


covariance function may be *rewritt en 


i 





jw(e_-e_). jwkT 
1 
t & Tae 
a =e 
y Mi oO 


E{y V8 


Now that the meanof the ensemble has been taken, one may apply the 


Mmit wand 


cof | pa 
Ki = | J ixyol® >) 


jw(e_-e€ 
E1y v9 


K(T) is now in the form of an inverse Fourier transform; thus the 
Wiener-Khintchine Theorem for two sided density functions may be uti- 
lized at this point, Let S'(f) be the two sided power spectral density. 
; rea ; 

jw(e jwkT 


) 
sine ixyol?) efyye 2 poe 8 
O 


k=-00 kT) 


-co < {f< © 


Later, the power spectral density shall be defined for positive fre- 
quencies only (one sided), and since S'(f) is an even function, the one 
sided power spectral density is S(f) = [ 2S'(f)] u_(f), where u (f) is 


the unit step function, The unit step function is defined as 


0 , £<0 
Umi ee eS 0 

S 
1 , £270 


Z0 





Define Sty (f) as the two sided spectral density for the typical sampling 


pulse, or 


Siz ff) is a continuous even function of frequency. The one sided spec- 


tral density for the typical sampling pulse would then be 


2 
2|X,,(£) | 
= CC 2 
S,,(f) = 25,,'(£) _ pe .0 
O 
= 0 iO 


Thus, the power spectral density may be rewritten 


im(e_ = iwkT 
jw(é,-e,) jWkT . 


CO 
! = ! = sfs 
SG) Su (f) ) E{yjy5¢ 7 , “co =f Fm. 


k=-00 ki 
O 


The spectral model to be used assumes that the y, 2re indepen- 
dent from sample to sample, and independent from the € process, 


Thus, the equation that is used for the generalized spectral model is 


nen 16 @) 

s(t) = Sf] y7 +7" +). 
k=-00 

10 


1) jokT | 
} e , “Co =f Soo, 


Observe that the model can also be used for pulse-amplitude modulation 
(PAM) in which the Yr represent the voltage level of the ae sample of 


the message and the pulse is subjected to a random jitter (in time), 


76:3 








possibly due to equipment temperature fluctuations or varying transmis- 
sion path length. The ae would of necessity be independent random 
variables for the generalized spectral model just developed, The eh 
however, as will be shown in the further development of this PPM 
model, can either be independent, or correlated random variables with 


a Gaussian distribution, 


(age 








De A Discussion of the Statistics in the Spectral Model 

As seen in the development of the spectral model, certain sta- 
tistics of the two processes € and y, together with the pulse shape, in 
essence determine the power spectral density. Since independence was 
assumed, the statistics of the pulse amplitude y, are quite simply de- 
termined, as all that is involved is the mean of the square of the pro- 
cess eo and the square of the mean Ge. The statistics of the € pro- 
cess are not as simple, especially since independence of ey will not be 
assumed for all cases. The statistics of the € process are directly re~- 
lated to the message process x(t). It has already been shown that the 
variance and the standard deviation of the two processes are linearly 
related. Furthermore, the probability density functions (pdf) of the two 
processes are of the same form if one considers continuous distribu- 


tions only (i,e., if x is Gaussian, € is also Gaussian). Specifically, 


P then 


2) 0 € 
p Ss fea 
Ee, (€) ve; PS (2 ; 
€ k eC 


Wo show that this 16 tries recall that c = (F /9 )x . Now 
x 


att P. (x) is the pdf for x 
k 


O € 
x 
< = < —— 
P[ eo e | Ie x 
€ 
But 
Ome O e/O 
P < —— = = 
x = P. (x )dx 
€ -0O k 


Since the pdf is the derivative of the distribution function, then 


Ze 





from which one can obtain the desired result 


oO O € 
‘ 8) x ‘e} 
k € k e 


Pe 


To avoid carrying unnecessary terms in equations already complex with 
terminology, atime axis shift of the PPM waveform will be assumed 
such that the mean of the process € is zero. That this simplification 
may be performed with no effect on the resultant power spectra is in- 
tuitively obvious, but nevertheless has been previously verified by the 
writer in each case considered. 

To evaluate the statistics of the process € that appear in the 
equation for the power spectral density, one must determine the joint 


characteristic function G(jw, jw) of a bivariate process, The joint 


Z 
characteristic function is defined as 
j(w, Ss jw E (e =) 
Glu nae Efe Z 2) ..: 2 1), 
J he J >) ~ j ° 
This joint characteristic function will be evaluated for the Gaussian 
bivariate case. For the remainder of the cases included in the compu- 
ter model, statistical independence from interval to interval of the ¢ 
process was assumed, For the independent case, the problem reduces 


to the determination of the square of the absolute value of the charac- 


teristic function of €, 
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me Efe} |S = la suwl” 


ef. 


The equation for the power spectral density now reduces to 


= CO 
wk T 
s(f) =Seyfy? + 77 |G w)|? ) eM - wo $f 500 
k=-00 


pit 
~ jwk T 1 = k 
JWK1g _ + )) ae pe 
de Tr sf Ts ) 
-0O G60 O 


Recall that Si, (f) = 2S'(f) for f #0 and S(f) = 2S'(f)u_(f). Then 


ae oO 
we? 2 il k 
oii) = S,,(f) y +y IG _(w)| -l) + 7 ) s( a ) u (£) 
oF k=0 O 


for the case in which c and y, are independent from sampling interval 
to sampling interval. As mentioned previously, only positive frequen- 


cies are considered, 


Zo 








5A, The Message as a Gaussian Process 

The purpose of this subsection (and the following subsections) is 
to fit a specific process into the generalized spectral model. The pro- 
cess considered may either be the message process x or the € cece 
as the similarity of the two processes has been previously shown. For 
the Gaussian case, the development will allow for a statistical correla~ 
tion to exist between I adjacent samples of the process, More speci- 
fically, this development will allow for the spectral model to handle 
Gaussian een eaSEER for which the autocorrelation function R(T) exists 
for x(t) up to |T| = IT |, and R(T) = 0 for | T | =iT 4 This ds practi- 
cal because a purely random process yields an autocorrelation function 
that approaches zero as T increases [1]. 

The expected value 


gfe 2} 


kT 


is the joint characteristic function of ¢ for a given kT =t.,-t 
O 


Dl 20. ae 
The joint characteristic function G(jw, » jw) aS previously defined fits 


the case now under development when w, = -w It was decided to 


] 2° 


carry out the general development of G(jw5 > JW), substituting into the 


final result w =w_ = #w 


3 i° Tinas- 


jue, - j(w, e+ é_) 
E{e 


eC 


The expected value is determined as 


pe G(jw), ju.) 
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2 ry 
G(jw, » jw.) -Jo fre p(e,,¢€,)de,de, 


The bivariate Gaussian pdf is 


[-1/20°(1-p')] (ec, “-2pe Ete >; 


| i }-. loa. 2 
ple),€,) = _ a >. 5 e 
270 V1 - p 
where 
Uo = O and O = an = BO 
] Z 
For simplicity, substitute 
ae: ca 
Z 
ZO MS l - p = b 


fiiaen 


1 00 jw ty /ab*) 
cumin) = ap fs 


= OO 
J) nN 


OO 
{ _-(/ab* )(r “- 2619). ava 


~ OO 


But 


CO Z jr 
| Bec eD ie -2ph¥) ad 


-0O 
lied oo fab™ fap] PU Gs 
-00 


which is in the form of a constant multiplying the characteristic function 


= pt and 0, © = ae 


of a single Gaussian process }\, with mean H. 


The characteristic function of a single Gaussian process may easily be 
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determined With the aid of tables of integrals), and the last integral is 
simply 


pa 
2.2, ,2f joyw,-(ab w, /4) 
i i 
Bbvat er y /ab e 


The joint characteristic function then becomes 


2 2 2 oN | 
lL -(abow, /4) ims -Wab™)y a ae 


GOO) ~ mea 


-OO 


Define 





and this integral also takes the form of the single characteristic function 


of the process yj, with the mean u, = 0 


y 
-. ; -(ab“w V/A ep (1/20, “0° jw 4 
G(jw) > jw) = syle e a, Ver \_ e & dy 


M2 ~(ab“w,“/4) (0, °w,°/2) 
O aj— e iS 
Vi a 


By resubstituting the values for Fy a, b, and W., one obtains 


1s Z Z 
=(6 /2)\(w, +2pw wo tw ) 
1s 2 


Recall that for this development, w = w, = -w 


3 i° Thus, 


Lee 
moe w (l-p) 
G(jw 1 jw) = G(w) =e 


Applying the above to the spectral model for a given value of kT: we 
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have 


22 
= -O l-p(kT 
c, 5 [l-p(kT )] 
=e 
kT =t_-t 
o 2 | 
That the expression is a function of the autocorrelation of the process ¢€ 
(or x) may be determined by observing that p is the correlation coeffi- 


cient for a given kT in the development. Considering all kT, p({kT ) 


is the normalized autocorrelation function, 


Now let R(T) = 0 for k >I. The power spectral density for this 


case is then 


2 
> > I -o ‘w [l-p(k)] jwkT 
Sf) = SMA, yo + ¥ )) 2 : 
k=-[ 
k#0 
(dae are 
-o Swe EEN) jw T © — jwkT 
+ e ( ) e te e ) 
k=-00 k=(I+1) 


But p(kT ) = p(-kT), as R(T) is an even function, and 


Z 
= ; I _o Sale atae eer 
S'(f) = Sy (4) yo ot 2) 5° “e , cos wkT | 
k=1 


Vater 
-O I = jwkT 
we -~ 1+ pe ce : 
k=-00 


and the spectral model for the Gaussian case (for f = 0) becomes 


"> : aoe “w°[1- p(k)] -o “w* 
S(f) = 5, (f) < 1 Te 2 DAs -e icos wkT 
: Zoe 
-O Ww pats 
aceeres © aro) | 


° K=0 
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aley, The Message as a Uniform Process 

The specific spectral model will be developed for the case in 
which the message determining the modulation of the PPM waveform is 
a uniform process, Statistical independence of both the € andthe y 
processes will be assumed, For this assumption, the spectral density 
was shown to be 

—_ | ec 
S(f) = S.(f){y° + 7° Ja tu |? | k = )) aC - oa bu (6) 
Oo k=f O 

The problem is therefore reduced to that of determining the character- 


istic function G(w) of a uniform process, 


CO ne 
G Z 
2) | e P, (e)de 


As shown previously, if x is uniformly distributed, then ¢ is also uni- 
formly distributed. For the case of uniform distribution of x, the total 
possible deviation e(max) allowed is TO/L seconds (or +(T /2L) from 
the zero mean assumed), Therefore, 


_ To To 

_ = ee Sec e peat 

P, \¢) See A. oe 
k Oo 


Thus, 


32s | 
Cae i O jwe 
-(T/2L) 


ot 


wT 
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sin Ty, 
wT 
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2L 


G (w) 
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The power spectral density for the uniform case then becomes 





wT é 
— sin coed 
2 Z Ze k 
= 7 _—_— al — ) § -— 
S(f) = S(f)¢y" + 9 ar + (: = 
Oo Ok) O 
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DG. The Message as an Exponential Process 

The specific spectral model will be developed for the case in 
which the message determining the modulation of the PPM waveform is 
an exponential process. Of interest for this case is that the speech 
process is approximately exponentially distributed, especially in the 
region of medium to large amplitudes [1]. 

For the exponential case, statistical independence will again be 
assumed for both the € and y processes. As was the case then in sec- 
tion 5B, the problem is reduced to determining the characteristic 
function of the process €, 


The exponential probability density function is 


1 -W2/o_)le| 


p_ (€) = = 
ae 4/2 oa 





from which one can determine 


= E = 0 
- (€) 


Vane) = ne) - UL = WB(e") = 0 


Define 


Then 
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The characteristic function is 


Om. 
G (w) = | el P. (e)de 
-0O k 


an (ee eure) © e(jw-a) 
5 | AS Je de + | Be a de 
-0O0 0 


Or 


2 
G @) = = 


S 2 2 
0) Gece 
The power spectral density of the PPM waveform for the exponential 
case is then 


oe 4 
2.82 6 k 
S(f) = Sylf y+ F ed - | Tr Y.5(1 : a u_(f) 
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Ge The Computer Model 

The PPM spectral model was developed in section four, and the 
model was adapted to specific modulation processes in section five, 

The resultant equations are quite formidable if manual calculation is re- 
quired to determine the power spectral density for various cases. For- 
tunately, the tedium of manual calculations of such proportions is not 
necessary nowadays, as the digital computer handles such tasks with 
great facility. A computer model will now be developed that will incor- 
po rate all cases previously discussed, plus the unmodulated case. 

The computer model developed relates the pulse positioning pro- 
cess ¢€ in terms of the original*modulating process, This allows the 
user to have a more direct insight as to the effects of the modulating 
process On the power spectral density of the PPM waveform. Mostof 
the parameters of the model are variable at the discretion of the user. 
An exception is the number of points of continuous spectra calculated 
between each possible discrete frequency spectrum. Since five points 
in each such interval are adequate for a graphic plot, this quantity has 
been fixed at five to avoid possible confusion by future users of the 
model. 

It should be noted here that the following development of the com- 
puter model may not at times appear to be the simplest approach. The 
approach taken generally yields simplification of computer Fortran oper- 
ations to be performed, such as eliminating the calling of a library 


function for performing an exponentation when a simple multiplying 
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operation wiil suffice. It is admitted, however, that a few of the ''sim- 
plifying" techniques employed did not necessarily make the program 
more efficient. The development will include these also, as they have 
become a part of the final computer model, 

While no extreme efforts were applied toward making the com- 
puter program efficient, this aspect was certainly kept in mind while 
the program was being developed, 

The program was written in Fortran IV language for the IBM 
7040 at Sylvania's Electronic Defense Laboratories (EDL). With some 
obvious card changes (for each library function used), the program can 
be adapted to the Fortran language of the CDC 1604 at the Naval Post- 
graduate School, 

The total run time (including card reading and program printout) 
for computing four separate cases in one program run from a binary 
deck was timed at two minutes on the 7040, 

From the previous developments, it was decided to base the 
computer model on the mathematical model expressed in the following 
form for frequencies greater than zero. The case for d-c (zero fre- 


quency) will be developed separately. 
= 2 
S(f) = Sylf){y-+ F[Sum(s) - |G(w)|7]} 


+ S (47° 1G (w) |e Yo(¢ d =) } 
Oe l O 


The first half of the expression is the continuous spectrum and the 


bes, 





second half is the spectrum concentrated at discrete frequencies. The 


term Sum/(f)} is defined as 


L, ~(o,0)"[1 pt] 


“(5 w) 
Sum(f) = #2 (. -e cos kwT : 
k=) 


The term Sum/(f) is included in the model so that the Gaussian correla- 
ted process may be handled, Note that when Sum(f) equals zero, the 
model reduces to that derived for the case in which all processes were 
assumed independent. The term G (w) is simply the characteristic 
function of the independent random process eé, and fits as shown into 
the chosen representation even when the Gaussian correlated case is 
considered. That this is true may be determined from the fact that the 
characteristic function for the independent Gaussian case (mean zero) 
is the eee root of that determined for the Gaussian bivariate case 


when rho (p) is zero, or 


SCONES 


(for Gaussian process with p = 0), 

All terms in the model have now been adequately defined except 
for the power spectral density of the "typical sampling pulse". Recall 
in a previous development that this pulse could be of any shape witha 
maximum amplitude of one. Also recall that this amplitude of one does 
not restrict the amplitude of the pulses in the PPM waveform, as the 
PPM pulse amplitudes are determined by the y process, Thus, the 


‘typical sampling pulse" is simply a description of the shape of the 
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pulses of the PPM pulse train. For this computer model, a rectangular 


pulse will be assumed, the pulse being of duration u and of unit ampli- 


tude, 


oud 
—_J,L 


Figure 6-1, The Typical Sampling Pulse for the Computer Model. 


As previously shown, 


2 2 
Sef) = IX, ,(£) | : 
O 
Here 
OO cate ‘euyi sae 
K(f) = | Ul(t)e 2 at = ye Sree 
~OO ~(u/2) 


from which one may determine 
Z 
eee 
ae aa = 
T Wu ‘ 
> pela 
Z 


5, (f) = 


To simplify the computer program (so that some of the time consuming 


operations are minimized), the above equation may be rewritten as 


2 
_ 4u— (1 = cos wu) 
Siy(f) eS 


os (gan 


The discrete frequencies at which S(f) is to be computed are 





To calculate S(0), one must apply the value of the unit step 
function {one~half) at zero frequency (the d-c term), Substituting w = 0 
into the equations developed thus far, one obtains for all cases 


S (0) 


Ue cene 2 eee 1 
s(o) = S— fy. 7?(1- 2) |} 
O 








where 
0 
a ue 
ie eG 
O 


The sequence of computation of the power spectral density may 
be determined from the flow chart (Figure 6-4). 

The computer model will determine two hundred values of S(f) 
up to and including f = 40/T. This includes all 40 possible values of 
S(f) concentrated at discrete frequencies (delta functions) within the fre~ 
quency range considered. 

A list of the significant symbols used in the program and their 
relationship to the terms in the mathematical model is shown in Table 
6-1. The computer program itself with one page of output is included 


as Figures 6-2 and 6-3. 
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Tables =i. 


Symbol 


A 

ALFA2 
CFSQ(N) 
CONT (N) 
DELEN 
1 


F(N) 


NR 


RHO(K) 
SIGMA 
S(N) 
SU(N) 


SUM(N) 


TEMP(K) 


bGsrel@) 


YooOr 


Computer Program Symbol Identification, 


Equivalent Symbol in Mathematical Development or 
Symbol Definition 


A 
a 3 
a (used for exponential case) 
2 
G (Ww 
}G_(w) | 
Value of continuous spectrum at f = N/5T 
Value of discrete spectrum component 
L 
Frequency; F(N) = N/5T. cps 
i] (upper summation index for Gaussian case) 


None; this index specifies modulation process, See 
comments in Figure 6-2, 


Number of separate cases to be computed this program 
run 


o{k) 


@) 
x 


S(f); Power spectral density 
S 
ff 
Summation for Gaussian correlated case 
Ak 
O 


Partial value of summation of SUM(N) 


ae a q 


Shy) 





PROGRAM PPM 
THIS PROGRAM CALCULATES THE POWER SPECTRAL DENSITY OF A PPM 
WAVEFORM THAT HAS AS ITS MODULATING SOURCE ONE OF THE FOLLOWING 
EITHER AN INDEPENDENT GAUSSIAN PROCESS OR A GAUSSIAN PROCESS 
WITH UP TO TEN CORRELATION VALUES(FOR TAU=T,2Tyeees10T)» AN 


~ JNOEPENDENT UNIFORM PROCESS, OR AN INDEPENDENT EXPONENTIAL PROCESS 
J=1 FOR GAUSSIAN,=2 FOR UNIFORM,=3(UNMODULATED) ,=4 FOR EXPONENTIAL 


AAA 


| 


ov. OV! 


C FORMAT 8 STATEMENT ALLOWS FOR READING RECORD IN FOLLOWING ORDER 
C U,sT,YBSQ, YSQB,SIGMA,A,FL,J:1,AND 10 VALUES OF RHO 


a ee ee — ee ee ee ei ee ee ee ee ee A a i a eae ae ae ere Te oe ~_- = — 


— DIMENSION RHO(10) ,SUM(200) »CFSQ(200),SU(200),5(200),TEMP(10) , 
1F(200), CONT(200) 
~~ FORMAT (4H F= FI0.3, 6H SUFI= EL2.4, SH CONTIF I= El2.4) 
7 FORMAT (/5X,1LHU,11X_1HTs 9X, SHYBSQ » BX 4HYSQB,y 6X» SHSI GMA, OX, LHAs10X, 
~ -L2HFL,4X,29H 2NO ROW IS Je, RHO(K=1,10)/) 
8 FORMAT (7F11.6/(213,10F6.3)) 
9 FORMAT (4H F= F10.3, 6H S(F)= E12.4) | 
PI= 3.14159265 
NR = 4 
300 PRINT 7 


| 
| 
| 
| 
I 
| 
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301 SUM(IN) = 0.0 
rr GOTO (,2939) 5S 8 
1 C= (0.4#PIT@SIGMA/(FL#A) )##2 
700 11 Net,200 === ~~ 
FLO2 = N 
— €C=C#FLO2eFLO2 
IF (CC.LT.23.0) GO TO 12 
CFSQUN) = 0,0 - 7 
GO TO 11 
12 CFSQ(N) = L.O/(EXPI(CC)) © 
11 CONTINUE 
TF (£.&Q. 0) GO TO 100 ~~ 
DO 18 N=1,200 
—  FLos = N° 
TEM = 0.0 
OO 15 K=1,!I 
FLO4 = K 
CD = C#FLO3*#FLO3*#(1.0-RHO(K)) 
IF (CO.LT.23.0) GO TO 14 
TEMP(K) = 0.0 ~ . - 
GO TO 15 
14 TEMP(K) = (L.O/(EXP(CD) )-CFSQUN) Je (COS10.4ePILeFLOZeFLO4) ) 
15 TEM = TEM*TEMP(K) 
18 SUM(N) = 2.08TEM 
GO TO 100 
2 PU = 0.28PI/FL 7 os 
00 25 N=1,200 
FLO6=N 
PUN=PU#FLO6 
25 CFSQ(N)=((SIN(PUN) )/PUN) ## 2 
GO TO 100 
3 00 35 N=1,200 
35 CFSQ(IN) = 1.0 
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Figure 6-2. Source Program for PROGRAM PPM, Page 1 
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GO TO 100 
4 ALFA2 = 2.0#(FL#A/(T#SIGMA) )##2 
PE = (0.4#PI/T)## 2 
DO 45 N=1,200 
FLO5 = N 
45 CFSQ(NI=(ALFA2/(PE®FLOS*#FLOS +ALFA2))## 2 
100 AA = 4.0#U#U/T 
= 0.4#PI#U/T 
DF= 0.2/T 
DO 110 N=1,200 
FLO] =N 
F(N)=FLOL # DF 
— G = BeFLO1 
110 SU(N) = AA#(1.0-COS(G))/(G#G) 
— SUZERO = AA/4.0 
SZERO = SUZERO#(YSQB -YBSQ#(1.0-1.0/T)) 
FZERO= 0.0 
DO 120 N=1,200 
120 SIN) = SU(N) #(YSQB+YBSQ#(SUM(N)-CFSQ(IN))) 
DO 125 N=5e200,5 
CONT(N) = S(N) 
DELFN = SU(N) #YBSQ#CFSQ(N) #1.0/T 
125 S(N) = CONT(N)+OELFN 
PRINT 9, FZERO, SZERO 
~ 00 210 K=1,196,5 
K3=K+3 
~ DO 200 N=K,K3 





N=K+4 
210 PRINT 6,F(N),S(N),CONTI(N) 
~ NR = NR-L 
IF (NR.GT.0O) GO TO 300 
~C THIS RUN FOR 10 KC SAMPLING. U = O.LT 
_ STOP 
a END +~— a 


SR on ee 








Figure 6-2, Page 2 
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START 


NR = # oi waveforms considered in this run| 


Read waveform 
characteristics 








rint waveform 
characteristics 


Set SUM(N) = 0] 
| N= 1,200 
J=4 / Sk 
(Gaussian) 







ential 
pagementia ) Commarea 


_GO TO J 









je 

(uniform ) 

Compute 
CFSQ(N) 













Compute 
CFSQ(N) 
N_ = 1,200 







FSQ(N) = 1,0 
=.1,200 











Compute 
SUM(N) 
N = 1, 200 






[Compute SU(N) 
(L oNp eet ye 


Compute SU(O) 
Com oni litle S(O) 


Compute continuous spectrum, CONT(N- 
itt SSU 200 






Compute discrete spectrum, DELFN 
for gaa N/ 58 


Compute total power spectral density 
S(N) = CONT(N) plus DELFN | 


Print F(N), S(N): also print CONT(N) 
when S(N) contains discrete spectrum 





No 
im NR = 02 iNR = NR = 1 
Yes 
Figure 6-4, Flow Chart for PROGRAM PPM, 
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a _ = 2 2. a ~_ 


ie Computed Power Spectral Densities, Results 

This section contains the displays of the power spectral densities 
computed for certain combinations of pulse widths, modulation processes, 
and varying degrees of pulse deviation (or RMS deviation from the mean 
position), The legend accompanying each graph indicates the pertinent 
data regarding the process for which that power spectral density was 
computed. 

The first graph (Figure 7-4) is the power spectrum of an unmo- 
dulated pulse train, Observe that with no modulation of the pulse train, 
there is no continuous component in the power spectrum, ~ 

Three different no emeiheed autocorrelation functions were used 
for the Gaussian case, For the process which has very little correlation 
between samples, the approximation used was the linear autocorrelation 


function indicated by Figure 7-1l, 


p(T) 


0 Dale ’ 
O 
Figure 7-l, Limited Correlation Between Samples. 


An approximation to the autocorrelation function of a process with a high 
degree of correlation existing for several adjacent samples is shown in 
Pigines! 2. 


AA 





1iT 


Figure. 7-2, Samples Highly Correlated. 


An approximation to the autocorrelation existing for band-limited Gaussian 
noise (0-W cps) is shown in Figure 7-3. In Figure 7-3, the Gaussian 
noise is assumed to have been sampled at ten times the highest frequency 


present, or at lOW samples per second, 


p(T) 


Figure 7-3, Approximate Autocorrelation forOversampled 


Band-limited Gaussian Noise, 


Two examples are shown for which the pulse amplitudes are not 
held constant (Rayleigh distributed), in addition to being pulse position 
modulated. It may be observed that in these cases the continuous spec~ 
trum begins with an appreciable value at zero frequency, which is a 


characteristic of pulse amplitude modulated waveforms, 





To allow for an adequate display of the continuous spectrum, a 
scale facto: was chosen that does not permit the showing of the discrete 


power spectrum to scale, The values of the discrete power spectral 


lines are numerically indicated for each line that exceeds the scale of 


the graph, 
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S. Spec ral Analysis of Actual PPM Waveforms 

A compar'son of the results from the spectral model was made 
with the results determined experimentally. The objective of this aspect 
of the investigation was to determine if short time interval analysis 
would yield a reasonable estimate of the power spectral density of PPM 
waveforms. 

In order to generate PPM waveforms comparable to those asso- 
ciated with the spectral model, a pulse position modulator was designed 
and built. The circuit and performance characteristics of the modula- 
tor are included in Appendix A. The unmodulated pulse train generated 
has a pulse repetition rate of one kilocycle, and a pulse width of 195 
microseconds, 

The spectrum analyzer (periodogram calculator) used for the 
experimental measurements is Known by the trade name ''Spectran", 

The Spectran is one of the parallel filter bank types of analyzers des- 
cribed in section 3. It contains 480 parallel filters, and has a resolution 
of 25 cycles per second (cps), The analysis time is 18 milliseconds. 
The first two frequency bands of the Spectran were used for the PPM 
analysis, covering 20 to 3852 and 3852 to 7684 cps, respectively. The 
use of two bands accounts for what appears to be a slight discontinuity in 
the pictures to be presented later in this section. 

The Spectran was calibrated by inserting a sine wave input of an 
arbitrary frequency (two kilocycles) and recording the analyzer response 


for measured power inputs, The input power was considered the square 
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of the RMS voltage input, which is equivalent to the spectral model defi- 
nition of power into a one ohm resistor. The resulting calibration curve 
is included as Figure 8-1, Unfortunately, the calibration curve is only 
accurate for spectra at discrete frequencies (delta functions) as a conse- 
quence of having been based upon data taken for discrete frequency in- 
puts. The calibration curve will not yield an accurate value of the con- 
tinuous Spectrum at a speciiic frequency fo. Each of the 480 "spectral 
windows" of the Spectran has a bandwidth of about 10 cps, Therefore, 
when a continuous spectrum is measured, the measurement obtained at 
a frequency Ee is actually the proportionate power of all frequencies that 
are passed through the spectral window that passes the frequency fi. 
Further, the noise level of the measuring equipment prohibited accurate 
measurements of the continuous portion of the spectrum. Consequently, 
for the experimental results, a determination of the exact values of the 
continuous spectrum must be attempted by a different procedure. The 
general shape of the continuous spectrum was determined by the Spectran 
analysis, however. and is displayed in Figure 8-3, A discussion of the 
background related to the photographs in Figure 8-3 is considered neces- 
sary in order that they may be properly interpreted, As predicted by 
the computer model, the value of the power spectral density at any fre- 
quency of the continuous spectrum is several orders of magnitude less 
than that of the larger discrete components. In addition, the predicted 
average levels of the continuous spectrum is about the same as the noise 


level output of the measuring equipment (the Spectran), Thus, if one 
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visually subtracts the noise level and the effects of the discrete compo- 
nent as determined from the top picture in Figure 8-3, then a fair esti- 


mate of the shape of the continuous spectrum may be gleaned from the 


bottom picture. 
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Figure 8-1, Spectran Calibration Curve 


A plot of (V at a discrete frequency vs Spectran out~ 


‘e 
RMS 
put response, (Spectran input attenuator set at 28 db. ) 


Nall, 2.0 S30 4.0 


Spectran Response in Volts 
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In order to measure the value of the continuous power spectrum 

at discrete frequencies, a wave analyzer with a tunable 6 cycle bandwidth 
i | ; 
was employed. The PPM waveform, modulated by a random noise 
(assumed Gaussian), was used as the input to the wave analyzer. True 
RMS voltage measurements of the analyzer response were made at 100 
; cy . ; 

cycle increments. The block diagram of the experimental setup is 


shown in Figure 8-4, 


lke 
iEulse fram 


G. R. Type 1390A 
Noise Gen. 
20 kc Bandwidth 
















|IPPmM 
Modulator 









Narrow band True Riis 
Tunable fates i Voltmeter 


HP-302A Fluke 910-A 





Figure 8-4, Technique for Experimental Measurement of 


Continuous Spectrum. 


Assuming an ideal rectangular response of unit amplitude for the 


filter, the total power to the RMS voltmeter is 


Le 
Power = | S(f)df , 
f -3 


O 


where f is the center frequency of the filter. It shall be assumed that 


l 

The analyzer is the HP-302A, which has an equivalent noise 
bandwidth of 6 cycles with the center frequency adjustable throughout the 
audio frequency range, 


Z 
True RMS voltage measurements were made with a ''Hi-damped" 
Fluke Model 910-A voltmeter, which is a power sensitive instrument 
calibrated in volts RMS, 





there are no extreme changes in the power spectral density in any 6 cycle 
interval: then Sif) may be approximated by a constant S(t). The meas- 


ured power is then approximately 
Measured Power = 6S(f.) : 


S(f._) may now be determined from the square of the true RMS voltage 
O 


measurement, which is the measured power (to a one ohm resistor), 





A plot of the power spectral density measured in the manner just 
described is shown in Figure 8-5, The few inconsistent measurements 
recorded are attributed to this writer's inability to ''time average" a 
slightly fluctuating voltmeter needle by visual observations. : i OTEOm— 
parison purposes, the theoretical plot of power spectral density is 
included as the dashed line in Figure 8-5, 

The measurement of the discrete spectrum was a simpler matter, 
as Output levels well above the equipment noise levels were being meas- 
ured. As previously mentioned, the calibration chart applies to the 
measured results of the discrete spectra, from which may be determined 
a reasonably accurate measured power level. Unfortunately, the smooth- 


ing of the Spectran ''readout"' is not perfect, resulting in the occurrence 


For a discussion of the variance of the output of a square law 
device (the same type as the measurement fluctuations encountered), see 
page 256, Davenport, W. B.. and W. L. Root, An Introduction to the 
Theory of Random Signals and Noise, McGraw Hill, 1958. 
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of peaks and valleys when one sweeps across the spectrum with a con- 
stant level) ESrput skenal. Thus, even the discrete frequency power mmeace 
urements will have some equipment measuring error introduced. 

The discrete spectrum displays are shown in Figure 8-2, From 
Figure 8-2 and theoretical computations, an experimental versus theore- 
tical comparison was made and is included in Table 8-1. The PPM de-~ 
vice is that described in Appendix A. modulated with a wide band Gaussian 


noise source, 


\ 
| Discrete Power Spectral Density 


t ets ces 
| i 
Discrete | 


t Unmodulated Pulse Train 


% 


Frequency 


¥ : ° ° 
il Experimental | ; Experimental 
| ame P : Thieormercalt aes 
Ii(From Fig. 8-2) | ( 
1 


Modulated PPM 














Z 
O Ewer & + = 2,72 
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ie 
From Fig, 8-2)| heoretical 







No zero 
frequency 
response 


No zero 
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Table 38-1, Experimental and Theoretical Values of 


Discrete Power Spectrum, 
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Oe Discussé#on of the Results and Conclusions 
A reasonable question to be answered now is ''What specific sta- 
tistics of the modulation process contribute to the shape of the power 


' This is one question that is answered in this section. 


spectral density °' 
Other items discussed in this section are: 

1, The ''general shape'' of the power spectral density of PPM 

waveforms, 

2. The power spectral density as a function of the modulating 

process, 

3, The parameters that determine the relative amount of the 

total power that is contained in the discrete and continuous 
Spectra. 

4, The effect of correlation between pulses in the pulse train 

with respect to their relalive positions Imgmime, 

5, Some general comments on the experimental determination 

of power spectral density. 

The shape of the power spectral densities of PPM waveforms 
appears to be generally of the form (sin ie if one considers only the 
frequencies beyond the first lobe. The portion beyond the first lobe has 
been observed to be mostly cont:nuous spectrum for the modulation pro- 
cesses considered, The reason for this can be determined from the equa- 
tion upon which the computer model is based. The characteristic function 


it will be assumed that all pulse waveforms discussed in this 
section consist of rectangular pulses except when otherwise stated. 
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for all moduiat’on processes considered ''damps out'' toward zero value 
as the frequency ‘s increased. resulting in a power spectral density at 
higher frequencies that ‘s essentially the same as the spectrum S(t) of 
the typical sampling pulse. 

The frequency at which the first lobe is terminated (where the 
envelope of power spectral density tirst becomes zero) is the reciprocal 
of the pulse width. The shape of the first lobe appears to be quite depen- 
dent upon the modulating process. Figures 7-5, 7-8, and 7-14 are 
representative of the spectrum of PPM waveforms that differ in the modu- 
lating process: i.e., a uniform process, an exponential process, anda 
Gaussian process, respectively. The biggest difference in these power 
spectral densities is in the distribution of power in the discrete frequency 
components. : A close look at the Figures 7-5, 7-8, and 7-14 addition- 
ally reveals a distinct difference in the continuous spectrum within the 
first lobe, As was the case with the discrete spectrum, this difference 
may be attributed to the type of statistical process used to modulate the 
meM wavetorn. 

The question asked at the beginning of this section will now be 
answered, For the cases in which statistical independence between 
samples is assumed. the square of the absolute value of the character- 
istic function is the only statistic of the modulating process that influ- 
ences the power spectral density of the PPM waveform. For all cases 

The frequencies at which the discrete components exist are 
multiples of the reciprocal of the sampling rate, or the average pulse 
repetition rate, 
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considered i: this investigation, the absolute value of the characteristic 
function is simply the Fourier transform of the probability density 
function of a process with zero mean, The process with a non-zero 
mean yields a characteristic function that is equal in magnitude to the 
characteristic function of the process with a zero mean, but with a phase 
difference. Taking the absolute value removes this phase difference, 
thus yielding the same results as for the process with zero mean. 

For the case in which statistical independence between samples 
is not assumed, the power spectral density is influenced by the joint 
characteristic function of the modulating process. The effects of corre- 
lated samples will be discussed in greater detail in the latter part of this 
section, 

The relative amount of the total power in the continuous spectrum 
of a PPM waveform with pulses of fixed duration and amplitude appears 
to be determined by the degree of fluctuation of the pulse shift. This ob- 
servation may be confirmed by comparing figures based upon like modu- 
lating processes but with different standard deviations (or different total 
pulse shift permitted in the uniformly distributed case), Comparing 
Eagpures ¢-l1l and 7-14, fomwhich on is 0.905 and 0.1, respectively, one 
may observe that an additional 21.6 percent of the total power is included 
in the continuous spectrum for the greater fluctuating pulse shift. For 
the uniform case, the PPM fluctuations may be increased to the maximum 
level by permitting the maximum pulse excursion (delay) to be equal to 


the full sampling interval Tv. less the pulse width u. This corresponds 





tO) lee de 1s u) in the PPM models developed. When this is done, 
the discrete spectrum becomes quite small, and eventually disappears 
for very narrow pulse widths. The latter observation may be confirmed 
by direct calculations with the mathematical model; the former observa- 
tion may be confirmed from Figure 7-5. In Figure 7-5, the maximum 
pulse excursion within a sampling interval is permitted. The pulse width 
SOs LT andsthe apeprecate powedwin the discrete components 15 Reduced 
to less than one percent of the total power (excluding the d-c component). 

Another factor which influences the relative amount of total power 
contained in the continuous spectrum is the pulse width. For a given 
modulation process, as the pulse width is decreased, the percentage of 
the total power contained in the continuous spectrum increases, In con- 
trast to the considerations of the preceding paragraph, in which the total 
power of the PPM waveform is constant (independent of the fluctuations), 
if the pulse width is decreased the total power decreases unless a corre- 
sponding increase in pulse amplitude is made. A plot of the percent of 
total power in each component as a function of the pulse width is shown 
in Figure 9-1 for a single modulation process (the independent Gaussian 
case with a= 0.05), Data for this plot was obtained from the computer 
moedel. 

The effect of statistical correlation between sample values of the 
modulation process tends to center the continuous spectrum about the 
discrete components. It is interesting to note that the relative amount of 


the total power in the discrete components remains essentially the same 


79 





for the correlated and the uncorrelated cases considered, yet the shape 
of the continuous spectrum is radically changed (see Figures 7-11, 12, 
and 13, for example), However, in the limit as the correlation function 
of the modulation process approaches unity, the PPM waveform is essen- 
tially modulated by a d-c level, and the power spectral density is all at 
discrete frequencies, In fact, the spectrum becomes that of an unmodu- 
lated pulse train (Figure 7-4), 

The correlated case is a practical consideration for typical PPM 
pulse trains encountered. For most PPM systems, there is time divi- 
sion multiplexing of several channels of data pulses, which are further 
interleaved with synchronizing (sync) pulses. The sync pulses are trans- 
mitted at a periodic rate. The exact periodicity of the sync pulses re- 
sults in a contribution to the discrete components of the power spectral 
density. Further, the data pulses are normally statistically related to 
the sync pulses, The mean position of the data pulse is related to the 
sync pulse by some fixed interval for each process. The degree of cor- 
relation between data and sync pulses is set by the data pulse variance 
about its mean displacement. Finally, adjacent data pulses may be 
correlated, as was the case for the Gaussian correlated process of the 
simpler single channel PPM model used for this investigation. 

The verification of the theoretical results by taking experimental 
measurements was discussed in section 8 The measurement of the 
continuous spectrum of the power spectral density can be quite arduous. 


The difficulty is partly due to the fluctuations of the modulating process 
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(or the degree of modulation) varying in a probabilistic manner during 
the measuring interval, [his effect, of course, can be minimgzedie, 
averaging the measurements over a longer measurement interval. Fur- 
thermore, the noise introduced by complex measuring equipment can 
obscure the low-level signals encountered in the continuous spectrum, 
The slight discrepancy between the measured and the calculated 
power spectral densities displayed in Figure 8-5 may be largely attributed 
to the non-linearity of the pulse-position modulator for the larger modu- 
lating voltages. Consequently, the standard deviation of the pulse shift 
€ is not as large as that determined from the measured standard devia- 
tion of the modulating process, This results in a smaller continuous 
component and corresponding increases in the discrete components as 
compared to the calculated power spectral density, which was based upon 


a linear modulation scheme, 
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APPENDIX A 
The Pulse~Position Modulator 

A pulse-position modulator was designed and constructed so that 
experimental measurements could be made on the type of PPM waveform 
used for the mathematical model, The design specifications were to 
build a device that sampled a modulating process at a periodic rate and 
generated a pulse train, the pulses of which being delayed in time by an 
amount linearly related to the sampled value. To accomplish this, it 
was decided to build a sampler which had as its "hold" circuit an RC net- 
work with a time constant that resulted in an approximately linear decay 
during the interval under consideration. Thislinearly decaying voltage 


''eero modulation'' reference, When the two 


was then compared to a 
voltages were equal, a pulse was transmitted toa wionestaple multivi- 
brator, The monostable multivibrator was used to give a pulse output of 
the desired shape, and to allow easy variation of pulse width and pulse 
amplitude. See Figure A-2 for pictures taken of the unmodulated and the 
modulated pulse trains at the multivibrator output. 

The circuit diagram is shown in Figure A-l, The variable param- 
eters of the circuit are the pulse width (change RQ), the pulse height 
(change Nese and the sampling rate (change Ry and pulse rate of the 


pulse generator), The sampling rate is, of course, the average pulse 


repetition rate of the PPM waveform. 
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Pulse train modulated with wideband Gaussian noise, 


Figure A-2, Performance of the Pulse-Position Modulator. 
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